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Abstract. Spirals are common in Nature: the snail's shell and the ordering of 
seeds in the sunflower are amongst the most widely-known occurrences. While 
these are static, dynamic spirals can also be observed in excitable systems such 
as heart tissue, retina, certain chemical reactions, slime mold aggregates, flame 
fronts, etc. The images associated with these spirals are often breathtaking, 
but spirals have also been linked to cardiac arrhythmias, a potentially fatal 
heart ailment. 

In the literature, very specific models depending on the excitable system of 
interest are used to explain the observed behaviour of spirals (such as anchoring 
or drifting). Barkley [4] first noticed that the Euclidean symmetry of these 
models, and not the model itself, is responsible for the observed behaviour. 
But in experiments, the physical domain is never Euclidean. The heart, for 
instance, is finite, anisotropic and littered with inhomogeneities. To capture 
this loss of symmetry, LeBlanc and Wulff [34, 35] introduced forced Euclidean 
symmetry-breaking (FESB) in the analysis. 

To accurately model the physical situation, two basic types of symmetry- 
breaking perturbations are used: translational symmetry-breaking (TSB) and 
rotational symmetry-breaking (RSB) terms. In this paper, we provide an 
overview of currently know results about spiral wave dynamics under FESB. 



1. Introduction 

The spiral is an integral part of Nature: it can be seen in a snail's shell, in the 
layout of a sunflower's seeds and in the path of a falcon on a hunt, to name but 
a few. These particular instances are fixed in space, but spirals can also evolve 
in time: hurricanes and galaxies are common examples that come to mind. It is, 
however, rather arduous to conduct experiments on the latter physical objects, for 
obvious reasons. 

On a smaller scale, where experiments are easier to control, spirals in evolu- 
tion have also been observed in excitable media such as heart tissue, slime-mold 
aggregates, the retina or certain chemical reactions (such as the famed Bclousov- 
Zhabotinsky (BZ) reaction). In these systems, waves propagate by 'exciting' a 'cell', 
which in turn 'excites' some of its neighbours before falling into a 'refractory' or 
'unexcitable' state, followed by a 'resting' state, ready to be 'excited' should the 
wave come its way again. 

These systems give rise to beautiful images (as can be attested to in figure 1). 
While this in itself might yield enough interest to study them, there is also (at least) 
one serious reason to do so: spiral waves have been linked to cardiac arrhythmias, 
i.e. to disruptions of the heart's normal electrical cycle [54,55]. Most arrhythmias 
are harmless but if they are 're-entrant in nature and [...] occur [in the ventricles] 
because of the spatial distribution of cardiac tissue [30, p. 401]', they can seriously 
hamper the pumping mechanism of the heart and so lead to death. As a result, a 
full understanding of spiral wave dynamics in these media becomes imperative. 
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Figure 1. Spirals in excitable media. On the left, spirals in a solution 
of the Ginzburg-Landau equation [26], on the right, spirals in colonies 
of the slime- mold aggregate Dictyostelium discoideum [3] . 



2. Historical Perspective 

Numerous experiments and simulations have been performed with excitable me- 
dia, see for instance [4,5,7,17,36] (a selected bibliography can be found in appen- 
dix A). The various 'spiral' motions that are observed are classified according to 
their tip path, an arbitrary point on the wave front that is followed in time, as can 
be seen, for instance, in figure 2. Some of the standard possibilities arc shown in 
figure 3. 

In the literature (see [6,41,53] for instance), specific systems of partial differential 
equations (PDE) have been used to attempt to explain the observed phenomena: 
for instance, the FitzHugh-Nagumo equations 

itj = — iii'^ — v) + Am 

Vt = <^{u + (3 - -yv) 

of cardiology, where f3 and 7 are model parameters, u represents an electric 
potential and v a measure of permeability, or the Oregonator 

u, = l{u-u^-fv^^) + Au 

Vt = u — V + D^Av 

for the Bclousov-Zhabotinsky reaction, where /, and q (small) are the model 
parameters, Dy is a diffusion coefficient and u and v represent the concentrations 
of certain chemical reactants. 

These two systems of partial differential equations are instances of a general class 
of PDE, given by 

(2.1) utix, t) = f{u{x, t)) + DAu{x, t), 

where x £ R^, u : x R™ is bounded and uniformly continuous, D is an 

m X m diagonal matrix and / : R™ M™ is some sufficiently smooth function. 
General systems of the form (|2.1|) arc called reaction- diffusion system (RDS).^ 



Winfree provides a very complete survey of their use as models [53]. 
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Figure 2. On the left, an isolated spiral in a RDS; the excited wave 
front is shown in red. On the right, the corresponding tip path. 




Figure 3. Possible paths of the spiral wave tip. From left to right: 
linear drifting, outer epicycle motion, inner epicycle motion and rigid 
rotation [35]. 



Four types of solutions of 12.1|l are of particular interest in the context of spiral 
wave dynamics: 

Rotating waves (RW) arc rigidly rotating periodic solutions that are fixed 
in a co-rotating frame of reference {i.e. the frame rotates uniformly with the 
same frequency as the solution). In physical and numerical experiments, 
the tip path is circular. RW are sometimes called vortices or rotors in the 
literature [53, 55]. 

Traveling waves (TW) are linearly propagating solutions that are fixed in 
a co-translating frame of reference (i.e. the frame translates linearly and 
uniformly with the solution). In experiments, the tip path of such a solution 
is a line. Strictly speaking, TW (or retracting tip waves [6]) are not spiral 
waves as they do not have a rotating component. 

Modulated rotating waves (MRW) are two-frequency quasi-periodic so- 
lutions that are periodic in a co-rotating frame of reference that rotates 
uniformly with one of the frequencies of the solution. The tip path of such 
a solution is a closed epicycle when the ratio of the frequencies is rational; 
otherwise the tip path densely fills a ring over time, with an epicycle-like 
motion. In the literature, MRW are sometimes called meandering waves. 
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Figure 4. Bifurcation diagram of spiral wave dynamics in 12.211 [5]. 



Modulated traveling waves (MTW) are rotating solutions, superimposed 
with a linearly propagating motion, that are periodic in a co-translating 
frame of reference that travels uniformly with the linear component of the 
solution. The tip path of such a solution is a helix-shaped two-dimensional 
curve. 

As an example of a RDS in which these occur, consider Barkley's system: 
(2.2) ^ " ' 

Vt ^ u - V, 

where a, b and c; are system parameters with s small [5]. Figure 4 shows a bifurcation 
diagram of the spiral dynamics of H2.2|l for £ = -g^j • There are three regions of interest 
labeled N, RW and MRW. Note that this last region is divided in two sub-regions by a 
curve labeled MTW. 

In N, no wave propagation is observed; in RW, observed solutions are RW and 
in MRW, observed solutions are MRW, with petality^ determined by the side of the 
curve MTW on which the parameters fall. The intersection of this curve with the 
boundary of MRW is a point that deserves special consideration: in every one of its 
neighbourhoods, the three basic types of spiral behaviours can be seen. 

An a priori surprising feature of reaction-diffusion systems is that figure 4 is a 
generic bifurcation diagram: most experimental results are strikingly similar [6,53]; 
this suggests they arc in fact a consequence of excitable media and their geometry, 
and not of the particular models that aim to describe the dynamics [4,21,34,35].'^ 
This is where the dynamical system approach enters the picture. 

^The orientation of the spiral "petals" . 

''it should bo noted that reaction-diffusion systems are not the sole models of excitable media, 
nor were they the first: Wiener and Rosenblueth originally defined and modeled excitable media 
using cellular automata [52]. 
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3. Equivariant Vector Fields and Reaction-Diffusion Systems 

Let r be a group acting linearly on a vector space X . A function / : X — > X is 
T — equivariant if it commutes with the action of F, i.e. 

Equivariant vector fields (with compact F) have been studied by many authors: 
notable amongst them are Golubitsky and Schaeffer [23], Golubitsky, Stewart and 
Schaeffer [24] and Vanderbauwhede [50,51] . The main feature of these F— equivariant 
vector fields is that whenever x{t) is a solution of x = f{x), so is jxit), for all 7 G F. 

The special Euclidean group SE(2) = C-j-S^ is a non-compact subset of all the 
distance-preserving transformations of the plane, with multiplication defined by 

(3.1) {pi,fi) ■ {P2,f2) = (e*'^'p2 +Pi,^i + ^2), V(pi,v3i),(p2,</'2) e §E(2). 

It acts on the space of bounded uniformly continuous functions from to M™, 
which we will denote by BC„(M2,R™), according to 

(3.2) ij-v){x) = {ip,^)-v){x)=viR_^ix~p)), V(p,^) e§E(2), 

where Rg represents a rotation by angle 9 around the origin. Reaction-diffusion 
systems on BCti(R^,M™) are §E( 2)— equivariant under the action of H3.2|l . but that 
action is not smooth over BCti(R^, R™): the problem arises with rotations, as a 
small shift in 9 produces a large displacement at far distances [56]. However, there 
is a closed set BCe(R^M'") C BC„(R2^R"') over which is smooth [56]. 

3.1. Abstract DifTerential Equations. In order to determine BCe(R^ M™), Wulff 
uses the following RDS paradigm [54]. Consider 

(3.3) ut{x, t) ^ DAu{x, t) + f{u{x, t),^), 

where x £ R^^ ^ . ]g2 ^ ]j+ ^ ]^m^ > is a diagonal matrix, <r G R^, A is the 
Laplacian and / is C'=+2 for some < k < 00. U detD ^ 0, let Y = BC„(R2,R") 
be the Banach space of uniformly continuous, bounded functions from R^ to R™. 
Otherwise, as long as / satisfies some additional growth conditions, the choice 
Y = L^(R^,R'") can also be used, with shght variations (see [27,44] for details).'' 

The semi-linear differential equation on Y associated to is the abstract 

differential equation 

du 

(3.4) _ = _Au + i^(u,i,0, 
where F{u, t, <;) ~ /(""(-, t), <;) and 

(3.5) A==diag(-diA,...,-d™A). 

Solutions of (|3.4|l are in one-to-one correspondence with solutions of (|3.3|l . 

In the remainder of this section, we assume the reader is familiar with basic 
definitions and results from operator theory and abstract differential equations (see 
[2,27,56] for details and definitions). 

Proposition 3.1. [56] Let A be given by (j^^jl, Y = BC„(R2,R'"), a G (i,l) 

and Ai = idy — DA, where D is as in Then A is sectorial in Y°' and 

■Tp-A7°', jr-^V" are bounded on Y. 
0x1 1 ' 0x2 ^ 

"^However, physical considerations demand that spiral waves in an infinitely extended medium 
be located in BC„ (M^.M™) [56]. 
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3.2. Existence and Uniqueness of Solutions. That 1)3. 4|l has classical solutions 
is shown by the following theorem. 

Theorem 3.2. [27, 56] Let Y , A, Ai and a be as in proposition 1^3. U be a 
subset ofV^ X R X M.^^ and F be as in ( [^.^t , locally Lipschitz in its first variable 
and continuous in the remaining variables. Then, for any (ug, tQ,<;) £ U , \ has a 
unique classical C^^"^ — solution u[t]UQ, to, <;) on [to,ti\, where ti = ti{uo,to,<;) > to 
and < k < oo is the smoothness of the nonlinearity f in 



3.3. SE(2)— Equivariance of Solutions. WulfF then shows that these solutions 
have a §E ( 2 )—cqui variant structure. A smooth local semi-flow {^t}t>o on a Banach 
space X is a smooth family of operators satisfying <I>o = id^ and 

(3.6) ^t+s = o $g = $s o $f for all s, t > 0. 

If furthermore $tx — > x as t ^ 0+ for all x G X, then {^t}t>o is a C^ — semi- 
group on X. A smooth semi-group on A" is a C*^— semi-group for which the map 
"ifx ■ (0, oo) X defined by = ^t {x) is smooth for aU x e X. 

The infinitesimal generator Lt of a smooth semi-group {Tt}t>o on a Banach 
space X is 

(3.7) Ltx = lim - ($tx - x) , 

t^o+ t 

whenever the limit exists. 

The special Euclidean group SE(2) plays an important role in the theory of 
spiral waves. For now, we assume it is parameterized as §E(2) = §0(2)-t-M^, with 
multiplication given by 

(3.8) (i?i, 5*1) • (i?2, <5'2) = {R1R2, Si + i?25'2). 

The standard E>E{2)-action on BC„(R^M™) is 

(3.9) {R, S) {v{x, t)) = V {R~^{x - S), t) . 

Note that §E(2) is generated by the families {5^}, {S^} and {Rg}, where 
„i /u\ fi2 /^0\ in (cosO —sin 6* 

satisfy H3.6|l . Of these, only {5^}^>o and {5'^}i/>o are smooth semi-groups on 
BC„(R^, M™); their respective infinitesimal generators are 

Lsi = --^ — and is^ = — • 
0x1 0x2 

On the other hand, {Re} is not a smooth semi-group (see [56, lemma 2.14]); the 
action of SE(2) on BC„(M2,M™) is not even continuous. Thus, BC„(M2,M'") is not 
a suitable space over which to define H3.3)l . 

This obstacle is overcome as follows. The formal evaluation of l|3.7)l yields 

d d 

Lr = X2-^ Xl- . 

0x1 0x2 

Now let Y = BCe(M^ M™) be the topological closure of domain(LH) in BC„(K2, M™). 
The reaction-diffusion system (13.311 is well-posed on Y, and proposition (|3.1|l and 
theorem 13. 21 still hold after substituting Y for Y [56]. Furthermore, (|3.9() is contin- 
uous on Y and the following result holds. 
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Theorem 3.3. [18,49,56] The semi-flow $t,<; generated by (l.^./jjl commutes with 
the restricted action ^SIB o/§E(2) over BCe(M^ M"). 

In particular, the other results of her thesis hold as long as ^ is a smooth 
§E ( 2 )—equi variant semi-flow (not necessarily generated by a RDS) on some suitable 
Banach space. 

4. Barkley's Insight 

Barkley was the first to realize that the Euclidean symmetry discussed in the 
previous section (as opposed to the specifics of a given model) could explain most of 
the spiral wave dynamics observed in experiments and simulations, and succinctly 
presented in figure 4 [4,5] . His key observation rests on the fact that for any reaction- 
diffusion system, the linearization at a RW at the onset of a Hopf bifurcation (hence 
at the boundary of MRW) has five isolated leading eigenvalues on the imaginary axis: 
= (due to rotational symmetry). At = (due to translational symmetry) 
and Xb = ±*/3o (responsible for the Hopf bifurcation from RW to MRW or vice- 
versa). 

The pairs of complex conjugate eigenvalues At and A^ can be made to coin- 
cide by varying two or more system parameters. The corresponding (interesting) 
codimension-two point is then found to lie precisely at the intersection of MTW and 
the boundary of MRW. 

4.1. Linear Stability Analysis at a RW. In his ground-breaking paper [4], 
Barkley considers the Oregonator-like system 

dui 1/1 f U2 + b 

— — = Aui + -Ui{l - Ml) Ml 

ot e \ a 

du2 r . 

— — = dAu2 + Ml - M2, 

ot 

where a,b,e are parameters with s small and S £ [0,1] constant. In vector form, 
this system may be written as 

(4.1) —=SAu + f{u), 

where m = (mi,M2)^, S = diag(l,(5) and / contains the remaining terms. The 
boundary conditions drU = is taken on a circle of radius R > 0. With this set- 
up, (|4.1|l is SE(2)— equivariant under the action of (|3.9|l . on some suitable Banach 
function space, as R —> oo. 

To find RW solutions of (|4.1() . i.e. solutions for which {dt + ujdg)u = for some 
speed of rotation it suffices to solve the eigenvalue problem 

F{u) =^ 

"-^■^^ DF{u)u ^ Am, 

where F{u) = 6Au + ujdgu + f{u) and DF{u) = 6 A + ujde + Df{u). Any A solving 
()4.2(l corresponds to an eigenvalue of the linearization of (|4.1|l at the RW solution. 

Using fast and efficient numerical methods, Barkley shows that three of the five 
leading eigenvalues^ lie on the imaginary axis. Indeed, the rotational symmetry 
of (|4.1|l forces Aj^ = (with corresponding eigenmode uji = dgu, where u is the 



'^Eigenvalues with largest real part. 
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spiral solution of the first equation in H4.2(l ^: the translational symmetry of (|4.1|) 
imposes At = ±iti^ (with corresponding eigenmode ut = dxU ± idyU, where u is 
as above). Note that this holds in spite of the fact that the boundary condition 
breaks the Euclidean symmetry: for sufficiently large domain, the real part of At 
is numerically indistinguishable from zero.^ 

There is a last pair of complex conjugate leading eigenvalues Xb = a{a) ± i/3(a) 
that crosses the imaginary axis for some prescribed a = a*, leading to a Hopf 
bifurcation or 'spiral wave instability', in which MRW are observed. 

All the remaining eigenvalues have negative real part and so do not affect spiral 
dynamics. As a result, the five leading eigenvalues are isolated in the spectrum and 
so any pair {u,Lu^.at) that solves 1)4. 2|l is not part of a continuum of solutions with 
continuously varying shapes or speed of rotation. 

These results are in fact model- independent; as long as the reaction-diffusion 
equations governing the field u are §E ( 2 )—equi variant, the five leading eigenvalues 
will have the above properties. 

4.2. The Ad-Hoc Model. Based on these observation, Barkley [6] constructed 
an ad hoc 5— dimensional system of ordinary differential equations (ODE) which 
replicates the above resonant Hopf bifurcation: 
p = V 

(4.3) v = v [fi\v\^,w^) + iwh{\v\^,w^)] 

w = wgdup, ui^) 
where p, u £ C, w G M and 

/(^,C) = -J + "i^ + "2C-C', .9(^,0 =e-C-l, and h{^X)=lo 

for some 70 G K. The variable p represents the position of the spiral tip, while 1; is 
its linear velocity and jqw its instantaneous rotational frequency rate. This system 
has RW solutions that undergo a Hopf bifurcation to MRW solutions and it also 
has a codimension-two resonant Hopf point. Furthermore, it is cquivariant under 
the distance-preserving planar transformations generated by 

/p\ /e'>\ /p\ /p + a + i(3\ 

R^\v\^\ e'-^v and T^,^ i> = v 

\w J \ w J J V ^ / 

where Rj and Ta^f^ represent respectively a rotation by angle 7 around the origin 
and a translation by the vector a + if3. 

Note the absence of p in the right-hand side of (|4.3() as position plays no role in 
Euclidean systems. Figure 5 shows the bifurcation diagram of l|4.3() for 70 = 5.6. 
The similarities with figure 4 are readily apparent, in particular when it comes to 
the presence of a codimension-two organizing center around which the three types 
of spirals can be found. 

5. The Dynamical System Approach 

Then, in what has been hailed a "major mathematical work on spirals [21]", 
Wulff [56] rigourously proved that the resonant unbounded growth observed by 
many authors (such as [4,6]) does indeed occur near the codimension-two point. 



^Barkley provides some very strong estimates to that effect. 
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Figure 5. Bifurcation diagram of llOt [6] 



The following result, the center manifold reduction theorem of Sandstede. Scheel 
and Wulff [19,43-45], remains, in the author's opinion, both the most technical 
general results on spiral wave dynamics and its most fruitful ally in applications. 
It is an extension to non-compact symmetry group of Krupa's [31] center bundle 
construction for relative equilibria and periodic solutions. 

5.1. The Center Manifold Reduction Theorem (CMRT). The CMRT helps 
provide a rigorous link between spiral solutions of H3.3|l and Barkley's ad hoc 
model (1131). Set se(2) = so(2) x where so(2) is the Lie algebra of SO(2), 
consisting of the 2x2 anti-symmetric matrices on M. As a one-dimensional vector 
space, 

' • 1^ 
-1 0, 



Define exp^^, 



(2) 



SO (2) = Span 
so(2) ^ §0(2) by 



= Span{ J2}. 



exP5o(2) (PJ2) 



cos b sin b 
— sin b cos b 



Then, se(2) is the Lie algebra of SE(2), when endowed with commutator and ex- 
ponential maps as defined by (|5.HI and l|5.2|l below. Let I2 be the 2x2 identity 
matrix. Then 



: (rir2 - r2ri, riS2 - r2Si) 
= (exp,„(2)(fi),f~^(exp3„(2)(rt) -/2)s), 
'2 J =0,1,2 [19,44]. 



(5.1) [(ri,5i),(r2,S2)] = 

(5.2) exp((r,s)t) = 

for i G M, Tj G SO (2) and Sj G 

Now, consider a RDS of the form H3.3|l . For a fixed a relative equilibrium of 
1)3.3(1 is a solution u{x,t) whose time orbit, or semi-fiow orbit, is contained in its 
group orbit under SE(2). More precisely, it satisfies 

(5.3) u{x,t) = exp((r, s)i)u(a;, 0) 
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for some (r, s) £ se(2). The isotropy subgroup of such solutions is 

= {cr G SE(2) : au{x,t) = u{x,t)}. 

If = 'Zi, then u{x,t) is an armed spiral. 

According to the definitions of section [3 RW are relative equilibria of (|3.3|l 
with no translation component; after a change of coordinates bringing the center 
of rotation to the origin, (|5.3|l becomes 

(5.4) u^{x,t) = (exp5j,(2)(r,i),0)u,(a;,0) 

for some non-trivial G SO (2). 

Similarly, a relative periodic solution of 1)3. 3|l is a solution 

u(a;, t) ~ exp((r, s)t)w(x, t), 

for some (r, s) G se(2), where w is a non-constant T— periodic function in t. Ac- 
cording to the definitions of section 12 MRW are relative periodic solutions of (|3.3() 
with no translation component; after a change of coordinates bringing the center 
of rotation to the origin, (|5.3(l becomes 

(5.5) u*{x,t) = (exp^„(2)(^*^)jO)'"^(2;i^) 

for some non-trivial r* G so{2) and a non-constant T— periodic function wJ The 
isotropy subgroup of u* is defined as for RW, and the interpretation is identical. 

Standard center bundle results (such as those presented by Krupa [31]) cannot 
be applied to H3.3|) because §E(2) is not compact (a small rotation at the origin will 
produce a large displacement away from the origin). However, it can be shown that 
relative equilibria and relative periodic solutions are members of BCe(M^, M™), and 
so, as was seen in sectional that the action of §E(2) defined by 1)3. 9|1 is continuous 
on RW and MRW. 

The two hypotheses that follow allow the resolution of some technical difficulties 
within the CMRT. 

Hypothesis 1. For the parameter <;i, [resp. assume [resp. u*] is a I — armed 
(normally hyperbolic) RW [resp. MRW\ of (I5'.51l with ^ r^, as in < \5.4i [resp. with 
^ r* and w as in (|5.,5)) ]. If D is singular, assume further that u* [resp. u*] is 
k + 2— times uniformly continuously differentiable. 

Scheel [46] has shown that such rotating wave solutions can arise from Hopf 
bifurcations in a large class of planar reaction-diffusion equations. 

Hypothesis 2. Assume that {/i : |/i] > 1} is a spectral set for the linearization 
exPBo(2)(-''*)-0$i, ?.(«*) [resp. exp5„(2)(-''*)'D$i,?*(u*)] and that 

dim(range(P,)) = 3, 

[resp. dim(rangc(P*)) ~ 5], where [resp. P*] is the spectral projection associated 
to li* [resp. u* ] . 

That this second hypothesis can hold has been verified numerically by Barkley [4] 
(see sectional). It should be noted, however, that Scheel [46] has also shown that 
this hypothesis fails to hold for a large class of asymptotically Archimedean spiral 
waves. 

^While we focus mainly on RW and MRW, the theorems of this section can easily be adapted 
for TW and MTW. 
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The following result establishes the existence of an invariant center manifold that is 
contained in an SE( 2)— invariant neighbourhood of the group orbit of [resp. u*]. 

Theorem 5.1. ([44], theorem 4, p. 142) For any close enough to [resp. ^*], 
there exists an §E{2) — invariant, locally semi-flow-invariant manifold M^^. Both 

and the action o/ §E(2) on Af™ are C'^+^ and depend C^^"^ — smoothly on 
Furthermore, Af™ contains all solutions which stay close to the group orbit of u* 
[resp. u*] for all negative times. Finally, M™ is locally exponentially attracting. 
Furthermore, the manifold is diffeomorphic to the bundle = C x §^ [resp. 
= C X T]. 

5.2. The CMRT and its Applications. In [21], Golubitsky, LeBlanc and Mel- 
bourne describe the structure of the equations on M^" assuming the spiral waves 
have trivial isotropy subgroup. Generally, the essential dynamics for Hopf bifur- 
cation from armed spirals are analyzed via a 5— dimensional system of ODE on 
the center bundle §E(2) x C describes. For £ = I, the general system reduces to 
the center bundle equations 

p^e''^FP{q,q) 

(5.6) ip = F^{q,q) 

q = F'^{q,q), 

where p,q eC, ip eE>\ F'^(O) = uj.^t G R, F9(0) = and DF'^iO) = ioj^,, is purely 
imaginary.^ The Euclidean action on SE(2) x C is given by 

(5.7) ix,9)-ip,ip,q)^ie:'p + x,ip + d,q), V(2;, 0) G §E(2). 

The analysis of (j5.6(l . the titular dynamical system approach, allows the authors to 
recover the results of Barkley and Wulff concerning the Hopf bifurcation from a 
RW and resonant growth by considering a parameterized version of (|5.6|1 : a quick 
note on Bogdanov-Takens bifurcation from 1— armed spirals is also provided. If 
£ > 1, the structure of the 5 dimensional center bundle equations changes, but, as 
an ODE system, it retains §E(2)— equivariance under the action 

(5.8) {x,e)„,-{p,^,q)^(e''p + x,^ + £9,e"'''>q), V(a;, 0) e §E(2), 
for some fixed m G {0, . . . , [^/2]}.^ 

Theorem 5.2. ([21], section 5, p. 571) Lety = N{y,fi) be the center bundle equa- 
tions for an £~ armed spiral, parameterized by fi eR. There is a unique parameter 
value iiq at which the spiral undergoes a codimension-two bifurcation to resonant 
growth if and only if £ and m are coprime. 

In a subsequent paper [22], the authors again use the dynamical system approach 
to show that while SO (2)— symmetry alone may explain rotating waves, Euclidean 
symmetry is necessary in order to observe the unbounded growth of Barkley and 
Wulff, as well as to explain the full bifurcation diagram of section |21 

Still, some observed spiral behaviours (sec next section for a partial list) arc left 
unexplained by a careful analysis of the center bundle equations (|5.6|l : this obstacle 
is overcome through the introduction of forced Euclidean symmetry-breaking. 



^The frequencies cJrot and ojper in 15.61 play similar roles to the parameters uj and /3o in [4]. 
'^This action is consistent with 15. 7i . 



12 



PATRICK BOILY 



6. The Effects of Forced Euclidean Symmetry-Breaking 

Physical experiments can never be perfectly Euclidean, if only because of their 
finite nature. In the heart, this reality is obvious. Cardiac tissue is anisotropic (i.e. 
heart fibres have a preferred orientation and electrical conductivity is direction- 
dependent). Furthermore, tissue distribution is not uniform: there are zones of 
relatively high density that affect cardiac activity [30]. 

Similarly, introducing light pulses in a light-sensitive BZ reaction changes the 
geometry of the system. Moreover, the boundary cannot be ignored when the size 
of the spiral core is 'comparable' to that of the domain. 

Yet, the heart and the BZ reaction retain a partly Euclidean local structure. At 
distances 'far' from the inhomogeneities, is their effect truly felt? If the anisotropy 
ratio is such that the 'preferred' direction is only slightly so 'preferred', are spi- 
ral dynamics really affected? Can the Oregonator distinguish the boundary from 
infinity if it is 'very distant' from the spiral core? 

The Euclidean model alone cannot explain these events: clearly, any model hop- 
ing to do so should incorporate /orcerf Euclidean symmetry-breaking (FESB) in order 
to maintain the 'partly Euclidean structure' described above. The combination of 
Barkley's approach with FESB predicts, amongst other, the following (observed) 
spiral behaviours: 

Spiral anchoring appears when local inhomogeneities are present: spirals 
are attracted or repelled by a RW which rotates around the site of the in- 
homogeneity. This has been observed in cardiac tissue [17] and in numerical 
simulations of a modified Oregonator [38]. 

Epicyclic drifting can be observed when the sizes of the physical domain 
and of the spiral core are comparable: the latter is then attracted to the 
boundary of the domain and rotate around it in a meandering fashion. This 
has been observed in experiments and numerical simulations in a light- 
sensitive BZ reaction [57,59]."'^'' 

Quasi-periodic anchoring is witnessed in periodically-forced RDS. The re- 
sults are similar to spiral anchoring, with the attracting/repelling structure 
consisting of either two- or three- frequency quasi-periodic motion. These 
have been observed in a light-sensitive BZ reaction which is periodically hit 
by light pulses and the corresponding modified Oregonator model [13,25]. 

Discrete RW and MRW can be seen in systems that incorporate the no- 
tion of anisotropy {i.e. the system has a 'preferred direction'). In general, 
the tip path has discrete two-fold symmetry. These have been observed in 
numerical experiments on the bidomain model of cardiac clectrophysiol- 
ogy [41]. 

Phase-locking can also be seen in systems that incorporate anisotropy: the 
rotation and meander frequencies can lock and this motion can be super- 
imposed with a slow drift. This has been observed in the bidomain model 
as well [41]. 

These are illustrated in figure 6. In what follows, we present the results obtained 
through systematic breaking of the Euclidean symmetry. 



^'^The term boundary drifting is also used. 

^^They are then called entrainment and resonance attractors, respectively. 

^"^Generically, these waves cannot occur in SE(2)— equivariant reaction-diffusion systems [7]. 
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Figure 6. Non-standard motions of the spiral wave tip. Top row (from 
left to right): boundary drifting, spiral anchoring, and a three-frequency 
attractor [35]. Bottom row (from left to right): a discrete MRW with 
Z2— symmetry and phase-locking with small drift [41]. 



6.1. A Single TSB Term. Using the center manifold reduction theorem of [43,44], 
LeBlanc and Wulff [35] showed that translational symmetry-breaking (TSB) from 
Euclidean symmetry generically leads to anchoring or quasi-periodic anchoring and 
that TW, MTW, boundary drifting and quasi-periodic attractors could also occur. 
They did so by studying a general perturbed ODE system on the center bundle 
§E(2) X C, taking the form 

p = e'^ [F^{q,q) + eG^ {pe~'^ .pe'^q^l.e)] 

(6.1) ^ = F^{q,q) + £G'^(pe-'^pe^^g,g,e) 

g = F'(g,g) + £G«(pe-'^pe'^g,g,e) 

where e G M is small and the G— perturbations are bounded and uniformly contin- 
uous in p and q. When £ = 0, H6.1|l is §E ( 2 )—equi variant under the action of 1)5.7(1 . 
but for £ 0, the system is generally only §0(2)— equivariant: the translational 
symmetry of the model has been broken. The following results are proved (directly 
or in equivalent forms) in [35]. 

6.1.1. Relative Equilibria. In the case of normally hyperbolic (rotating) relative 
equilibria, ^'^ we may assume without loss of generality, and after an appropriate 
time-rescaling of the ip variable, that the center bundle equations H6.1|) take the 
form 

(6.2) p = e'* [v + eH{pe-'\pe'\e)\, 

where v G and £ G M is small. Set H{w,w,e) = H{w — iv,w + iv, s). 



'That is, q = in 16. li and the RW n« is not at the transition to a MRW. 
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Theorem 6.1. Let a = Kc[DiH{0,0,0)]. If a ^ 0, then for all e small enough, 
the center bundle equation H6'.ijjl has a unique smooth branch of periodic solutions 

(6.3) p,{t) ^ {-IV + 0{e)) (p{t) ^ t, 

whose stability is exactly determined by the sign of as. 

These periodic solutions are centered around the origin in the p-plane and are 
observable as anchored RW in the physical space. Note that the hypotheses of 
theorem 16. II are generic. 

Theorem 6.2. Let 



I{p) = Re 



/ e-''H {pe~'\pe'\0)dt 
Jo 



If Pa > is a hyperbolic solution of I{p) ~ 0, then for all e =/= small enough, 
the center bundle equation (|6'.ii|) has an epicyclic solution around the origin, whose 
stability is exactly determined by the sign of el' [po). 

These solutions represent quasi-periodic motion around the origin in the p— plane 
and are observable as epicycle-like motion along a circular boundary in the physi- 
cal space, with angular frequency given by 1 + 0(e). Note that the hypotheses of 
theorem 16. 21 are not generic. 

Since TW can also be seen as oo— centered orbits with no rotational component, 
an appropriate change of variables considered with the limiting case w^ot — > takes 
the center bundle equations H6.1() to the equivalent system 

(6.4) i = -vz^ + iezCf{z,z,e) - ez^C''{z,z, e), 

where C"^(2,2,e) = G'^iz~\z~\e) and CP{z,z,e) = Gp{z-\z-\s). 

Theorem 6.3. If and are sufficiently smooth near z ~ 0, the center bundle 
equation ^6.4^ undergoes a transcritical bifurcation of equilibria at e = 0. 

The trivial equilibria z = represent traveling waves in the p— plane and are 
observable as linear drifts in the physical space. 

6.1.2. Relative Periodic Solutions. In the case of normally hyperbolic (rotating) 
relative periodic solutions, we may assume without loss of generality that, after 
an appropriate time-rescaling and a change of variables, the center bundle equations 
(|6.1|) take the form 

w = —iUrotW + ell^{w,w, t, e) 

If = Wrot + eli^[w,w,t,e), 
where H'^ and H''^ are 27r— periodic in t, oj^ot ^ ^ and w = pe~'^'f . 
Theorem 6.4. Let h^{t) ^ D^H'"{0, 0, i, 0) and set 



/3 = Re 



27r 



h'^{t)dt 







If [3 ^ Q, then for all e ^ Q small enough, the time— 27: map of the center bundle 
equations 1] 6'. 5|) has a unique smooth branch of hyperbolic fixed points whose 
stability is exactly determined by the sign of ep. 



^^That is, the q equation in 16. li has a 27r— periodic solution and the MRW u* is not at the 
transition to a RW. 
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These fixed points represent periodic solutions centered around the origin in the 
w-plane and are observable as anchored MRW in the physical space. Note that the 
hypotheses of thcorem l6.4l arc generic. 

There are certain similarities between theorems 16.11 and 16.41 in the same vein, 
the following two results arc related to theorem 16. 21 



Theorem 6.5. If uj^ot ^ Q, 

J(p) = Re 



1 



If Pq > is a hyperbolic solution of J{p) ~ 0, then for all e ^ Q small enough, the 
equation in (|6'. 5\ has a unique smooth branch of hyperbolic invariant two-torii 
Tg, whose stability is exactly determined by the sign of eJ'{po). 

Such an invariant two-torus represents an 0{e) drift of a MRW (centered at 
a point different from the origin) around its SO (2)— orbit about in H6.5|l and 
is observable as a three-frequency motion in the physical space. Note that the 
hypotheses of theorem 16. 51 are not generic. 

Theorem 6.6. If LOrot G Q, with oj,.ot — j^, gcd{q,j*) = 1 and f > 1, let 

T^oo T Jq 

IfS,(t) is a hyperbolic periodic solution [resp. an equilibrium point ] of the Zj. — egui- 
variant ODE ^ = /io(CjC); then for all e ^ small enough, the w— equation in 
(|6'. 5|) has a unique smooth branch of hyperbolic invariant two-torii T{e) [resp. of 
hyperbolic "^—periodic solutions Weit)], whose stability as an invariant set is exactly 

determined by the product of the sign for the stability of £^{t) with the sign of e. 

The interpretation of these invariant two-torii and/or periodic solutions are ex- 
actly as in the remarks following theorem 16 . 41 and 16.51 

MTW are approached much the same way as TW were tackled in the preced- 
ing section. After an appropriate change of variables and a time-rescaling, and 
considering the limiting case dj^ot ^ 0, the relevant equations are 

z = -vz^ + eE'{z,z,t,e) 
<^^eE^{z,z,t,e), 

where E''(z,z,t,e) = -z'^H'" {z-^ .z''^ ,t,£), Ef{z,z,t,e) = m{z''\z-\t,e) are 
27r— periodic in t and w G . 

Theorem 6.7. If E'^ and E"^ are sufficiently smooth near 2 = 0, the time— 27: map 
of the z— equation in the center bundle equations (|6'. 611 undergoes a transcritical 
bifurcation of fixed points at e = 0. 



The trivial fixed points z = represent MTW in the plane and are observable 
as linear meandering in the physical space. 
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6.2. A Single RSB Term. LeBlanc [34] then showed that rotational symmetry- 
breaking (RSB) from EucHdean symmetry could provide an explanation for the 
appearance of discrete RW, discrete MRW and phase-locking in excitable media 
with anisotropy. He did so by studying a general perturbed ODE system on the 
center bundle SE(2) x C, taking the form 

p = e^^ [FP{q,q)+eG^{^,q,q,e)] 

(6.7) ^ = F'^{q,q) + eG'P[v,q,q,e) 

q^F''{q,q)+eG'^{^,q,q,e) 

where e G R is small, j* £ N, and the G— perturbations are bounded and uni- 
formly continuous in p and q, as well as ^—periodic in ip. When e = 0, (|6.7() is 
SE(2)— equivariant under the action of H5.7|l . but for e ^ 0, the system is generally 
only C-j-Zj. —equivariant: the rotational symmetry of the model has been broken. 
The following results are proved (directly or in equivalent forms) in [34]. 

6.2.1. Relative Equilibria. Let j* > 1 be a fixed integer. In the case of normally 
hyperbolic (rotating) relative equilibria, we may assume without loss of general- 
ity that, after an appropriate time-rescaling of the ip variable, the center bundle 
equations H6.7(l take the form 

ip = cjrot + sG'^{(p,e), 
where v € , oj^-ot G K and C, G"^ are —periodic in ip. 

Theorem 6.8. Assume 7^ 0. For all e ^ Q sufficiently small, the solutions of 
(|6'.ij^ are "^—periodic in time with discrete "Lj-— symmetry. 

These periodic solutions represent discrete RW in the physical space. 

Theorem 6.9. Assume uj^ot = 0. For all e ^ Q sufficiently small, if G''^ {<p,0) ^ 
for all G the solutions of H6'.6^ are discrete RW, with (large) radii of the order 
of -. On the other hand, if there exists p* £ [0, 27r) such that G^{(p*,Q) — and 
DipG^{ip* ,0) ^ 0, then (|6'.(j^ has at least j* stable (attracting) TW solutions and 
an equal number of unstable (repelling) TW solutions. 

In the latter case, all solutions of H6.8|l end up drifting linearly, after an initial 
transient period. 

6.2.2. Relative Periodic Solutions. In the case of normally hyperbolic (rotating) 
relative periodic solutions, footnoteWhere the corresponding 27r— periodic solution 
q*{t) to the q equation in l|6.7|l is such that F"^ {q* (t) , q* (t)) ^ for all t. we may 
assume without loss of generality that, after an appropriate time-rescaling and a 
change of variables, the center bundle equations (|6.7|l take the form 

p = e''^[FP{0)+eGP{ip,t, e)] 

(6.9) V : 

^^uj,,,+F'^{9)+eG'^{p,t,e), 

where all functions are 27r— periodic in 6 and ^—periodic in ip, and where the 

average value of F"^ is 0. In the unperturbed case, all solutions of (|6.9|l are MTW 
(if Wrot e Z) or discrete MRW (if u,,, i Z). 
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If e 7^ is sufficiently small, the cp— equation in (|6.9|l defines a Zj* — equivariant flow 
on a two-torus, with associated Poincare map 

(6.10) P{lp] oijot, = + 27rcj,ot + £H{lp, t^.ot, e), 

where H is ^—periodic in if and C" for some r > 3. The map P is thus a circle 
map; denote its rotation number by Pu^^^.e- 

Theorem 6.10. If Purot,!^ — ^^^^^ 7 > 0, m G Z are coprime, set k ~ gcd(7, j*). 
Then, phase-locking occurs in (|6'.i^ when (iUrot,^) lies in the m : 7 Arnol'd tongue 
of P. 

If fc ^ 1, the p— component of solutions of 1)6.9(1 is a 7— pctaled 2777— periodic 
curve with Z^^— spatial symmetry. Otherwise, the p— component of solutions of (|6.9|) 
is a superposition of a motion akin to the one in the case k ^ 1 together with a 
'slow' linear drift. 

Theorem 6.11. If Puj^^-^.e ^ Q, the flow on the two-torus described above is ergodic. 
If there exist cr G (0, 1) and K > such that 

k 



J 



> K\j\ 



-{2+cr) 



for all non-trivial integer pairs (fc, j) (which is almost always the case), then the 
p— component of solutions of ((6.9() is quasi-periodic and the closure of its positive 
image has Z^. —rotational symmetry. 

6.3. Simultaneous TSB Terms. The next logical step lies in studying the effects 
of n simultaneous TSB perturbations, for n > 1, which is done in [8,9,11]. Any 
excitable media which is littered with inhomogeneities, such as the human heart, 
could then in theory be modeled by a general system of ODE on the center bundle 
§E(2) X C, taking the form 

n 

n 

(6.11) f = F^{q,q) + Y,^^Gf{{p-C^e-'^,'{p^)e"',q,qA) 

n 

q = F'^{q,q) + ^ A,G?((p - 6)e-'^ (^^e^^ g, 9, A) 

i=l 

where A = (Ai,...,A„) e R" is small, Ci,...,^n G C are all distinct and the 
G— perturbations are bounded and uniformly continuous in p and q. 

When A = 0, (|6.11|) is §E(2)-equivariant under the action of H5.7|) : when A ^ 
is near the origin and along the j^^ coordinate axis of M", H6.11|l is generally only 
§0(2) -equivariant (i.e. it only commutes with rotations about the point ^j), and 
when two or more of the Xi are not zero, the system is generally only trivially 
equivariant: the translational symmetry of the model has been broken. 
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6.3.1. Relative Equilibria. In the case of normally hyperbolic (rotatmg) relative 
equilibria, we may assume without loss of generality that, after an appropriate 
time-rescaling of the (p variable, the center bundle equations H6.11|l take the form 



(6.12) 



+ > A,i^,((p-^,)e-^^(p-0K^A) 



where, ?j G and the functions Hj arc smooth and uniformly bounded in p. Boily, 
LeBlanc and Matsui study spiral anchoring in this particular setting [8, 11]. 

A 27r— periodic solution p\ of H6.12|l is called a perturbed rotating wave of H6.12|l . 
Define the average value 

1 f'^"" 
[pa] A ^ ^ J P^{t)dt. 

If the Floquet multipliers of px all lie within (resp. outside) the unit circle, we 
shall say that [pa]a is the anchoring (resp. repelling, or unstable anchoring) center 
of PA. 

Theorem 6.12. Let k e {!,..., n} and define ak — Kc[DiHk{iv, ~iv,0)]. If 
Uk ^ 0, there exists a wedge-shaped region of the form 

Wk = {(Ai, . . . , A„) G M" : < Wk,j\Xk\, Wk,j > 0, for j k and Afc near 0} 

such that for all ^ \ E Wk, the center bundle equation t^fl.lS^ has a unique 
perturbed rotating wave S^, with center [S^]a generically away from ^k, whose 
stability is uniquely determined by the sign of XkOtk- 

In contrast to theorem 16. II these periodic solutions are not necessarily centered 
around an inhomogeneity £^k in the p— plane, but they are still observable as an- 
chored RW in the physical space. Note that the hypotheses of theorem 16.121 are 
generic. Furthermore, [S^]x is a center of anchoring when Afca^ < and a center 
of repelling when XkO-k > 0. 

Boily then showed that theorem l6.2l has a similar generalization [8,9] . An epicycle 
manifold of (|6.12|) is an invariant set £\ for H6.12|) in which all solutions are epicycles 
when projected upon the p— plane. 



Theorem 6.13. Let k E {1, . . . , n} and 

/ e-'*Hk {pe-'*- iv, pe'*+ iv, O) dt 



Ik{p) = Re 



If Pk > is a hyperbolic solution of Ik{p) = 0, then there exists a wedge-shaped 
region of the form 

Vk = {(Ai, . . . , A„) e R" : |A,| < 14.j|Afc|, 14,, > 0, for j^k and Afc near 0} 

such that for all X G Vk , id-li^ has an epicycle manifold £^ whose stability is 
exactly determined by the sign of Xkl'kiPk)- 

These solutions represent quasi-periodic motion centered near (but not generi- 
cally at) the inhomogeneity ^fc in the p— plane and are observable as epicycle-like 
motion along a circular boundary in the physical space. Note that as was previously 
the case, the hypotheses of theorem 16.131 are not generic. 



SPIRAL WAVES AND THE DYNAMICAL SYSTEM APPROACH 



19 



6.4. Combined TSB and RSB Terms. Yet another way in which the Euclidean 
symmetry can be broken lies in the combination of TSB and RSB terms; such a 
situation is analyzed in [8,10]. Anisotropic media near an inhomogcneity. such as 
cardiac tissue in the neighbourhood of site of higher density at the origin, could then 
in theory be modeled by a general system of ODE on the center bundle §E(2) x C, 
taking the form 

p = e^'fi [FP{q,q) + eGP{ip,q,q,e,fi) + fiRP {pe-^'^ ^pe''^ , q,q, e, fi)] 

(6.13) ip - F^(g, q) + £G^(^, q, g, £, ^i) + Mi/^(pe-«^, pe''^, q, q, e, ^l) 

q = F'' (q, g) + £G« (</?, g, g, e, /i) + ^lH^ {pe-^^ , pe"^ ,q,q,e,ii), 

where (e, (x) g is small, f G N, the G, i/— perturbations are bounded and 
uniformly continuous in p and (/, and the G— perturbations are ^—periodic in Lp. 
Throughout this section, we fix f G N. 

Let C+Zj- be the subgroup of SE(2) containing all translations and rotations by 
angle k £ Z. When (e,/i) = 0, H6.13|l is §E(2)-equivariant under the action of 
()5.7|l : when e 7^ is small and /i = 0, it is C+Z^. — cquivariant; when fi ^ is small 
and e = 0, it is SO (2)— cquivariant, and it is generally only trivially cquivariant 
otherwise: both the translational symmetry and the rotational symmetry of the 
model has been broken. 

In the case of normally hyperbolic (rotating) relative equilibria, we may assume 
without loss of generality that, after an appropriate time-rescaling of the ip variable, 
the center bundle equations H6.11|l take the form 

(6.14) p = e** [v + eG{t, e, fi) + fiHipe''* ,pe'\e, /i)] , 

where (e, /x) G M^, w G C, G and H are smooth and uniformly bounded in p, and 
G is periodic in t. The following results are proved (directly or in equivalent 
forms) in [8, 10]; they depend greatly on the nature of j*. 

6.4.1. The Case j* = 1. As G is then 27r— periodic in t, it can be written as the 
uniformly convergent Fourier series 

(6.15) G(i,e,/i) = 5]5n(e,/^)e"*. 

Theorem 6.14. //g_i(0,0) 7^ and if a = Re[DiH{-iv,iv,0,0)] ^ 0, there 
exists a wedge-shaped region of the form 

W = {(e,Ai) e : |e| < K\fi\, K > 0, n near } 

such that for all (e, /i) G W with e ^ 0, the center bundle equation has a 

unique hyperbolic discrete rotating wave Vl ^ with trivial spatio-temporal symmetry, 
centered near but generically not at the origin, whose stability is exactly determined 
by the sign of fici . 

Theorem 6.15. Let 



R{p) = Re 







e"**iJ (pe^**- iv, pe**+ iv, 0, O) dt 



If Po > Q is a hyperbolic solution of R{p) ~ 0, then there exists a wedge-shaped 
region of the form 

V {(e, Ai) G : |e| < K\^i\, K > 0, fi near } 
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such that for all (£,/i) G V with e 0, the center bundle equation ]6.14i has an 
epicycle manifold Ql ^ centered near but generically not at the origin, whose stability 
is exactly determined by the sign of fiR'{pQ). 

To get the complete picture (in both of these theorems), the situation would 
also need to be analyzed near the e— axis. This implies dealing with fixed points of 
maps at oo; as of now, it has been relegated to a future investigation. Note that 
the hypotheses of theorem 16. 141 are generic, while those of 16. 151 are not. 

6.4.2. The Case j* > 1. In fully anisotropic media, one would have j* = 2. 

Theorem 6.16. Let ci be as in theorem \6.14\ If ci ^ 0, then there exists a small 
deleted neighbourhood of the origin, such that for all ^ (e, /x) €: , the 

center bundle equation ^6.141 has a unique hyperbolic discrete rotating wave 'D^ ^ 
with Wjj* — spatio-temporal symmetry, centered at the origin in the p— plane, whose 
stability is exactly determined by the sign of fici . 

Theorem 6.17. Let R{p) be as in theorem \6.15[ If po > is a hyperbolic solution 
of R{p) — 0, then there exists a small deleted neighbourhood of the origin, such 
that for all (e,/i) G , the center bundle equation ^6.141 has an epicycle manifold 
Ql ^ centered at the origin, whose stability is exactly determined by the sign of pci . 

Contrary to what might be thought at first, the epicycle manifolds of theorem 
l6.17l do not have Zj* —spatio-temporal symmetry; however, the epicycles themselves 
possess this symmetry in an appropriate frame of reference. As has been the case 
throughout this article, the hypotheses of theorem 16 . 1 61 are generic, while those of 
theorem 16 . 1 71 are not. 



7. Conjectures, Related Results and Future Work 

The application of FESB to the study of rotating waves has yielded a number 
of interesting verifiable results: the predictive power of the dynamical approach 
cannot be denied. Yet the full picture of spiral wave dynamics is far from complete 



7.1. Conjectures. In this section, we present some conjectures concerning modu- 
lated rotating waves. 

7.1.1. Simultaneous TSB Terms. The equations describing the essential dynamics 
of a normally hyperbolic modulated rotating wave are similar to those of rotating 
waves. Near such a MRW, the center bundle equations H6.7|l are equivalent to 



p ^ e^ 



(7.1) 



where A 



^A,ffJ((p-e,)e-'^(p-?,)e'^i,A) 



)e-^^,{p~ij)e^^,t,\), 



(Ai, . . . , A„) G R", f G C, Wrot 7^ and the functions H^''^ arc smooth, 
uniformly bounded in p for all j = 1, . . . , n and 27r— periodic in t. 

System (j7.1(l cannot be analyzed as easily as (|6.12() . but it seems nonetheless 
likely that the following hold. 
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Conjecture 7.1. Let k E {1, • ■ • ,n}. Generically, there is a wedge-shaped region 
of the form 

Wk = {(Ai, . . . , A„) e R" : |Aj| < WkM, Wk,j > 0, for j ^ k and Afe near } 

such that for all ^ X E Wfc, the center bundle equations ( |7. J| ) have a unique 
perturbed modulated rotating wave solution centered generically away from 

Conjecture 7.2. Let k € {1, . . . ,n}. Given a hyperbolic equilibrium of a related 
averaged system, there is a wedge-shaped region of the form 

Vk = {(Ai, . . . , A„) e M" : |A,| < 14.j|Afe|, 14,, > 0, for j^k and A^ near 0} 

such that for all ^ X E Vk, the center bundle equations ( |y.i| ) have a hyperbolic 
^--frequency epicycle manifold £^ centered near, but generically not at, ^fe. 

7.1.2. Combined RSB and TSB Terms. Fix f e N. The equations describing the 
essential dynamics of a normally hyperbolic modulated rotating wave are similar 
to those of rotating waves. Near such a MRW, the center bundle equations H6.7|) 
are equivalent to 

p = e»^ [i; + £GP(^,i,e,M)+M^^(pe-'^pe'^i,e,/i)] 

^^Lo,,,+eG^{^,t,e,ii)+^iH'fi{pe-''^,pe''^,t,e,iJi), 

where (e,/x) G R^, w G C, Wrot 7^ 0, the functions G'^^'^^HP'^ are smooth and uni- 
formly bounded in p and 27r— periodic in t, and the functions G^'*^ are additionally 
^—periodic in ip. 

System (|7.2|l cannot be analyzed as easily as (|6.14() . it seems nonetheless likely 
that the following hold. 

Conjecture 7.3. Let j* = 1. Generically, there is a wedge-shaped region of the 
form 

yy ^ {{e, /x) G R^ : |e| < A>|, K > 0, ^l near } 
such that for all (e, /i) G W with e ^ 0, the center bundle equations \7.'^ have a 
unique hyperbolic discrete modulated rotating wave ^ , with trivial spatio-temporal 
symmetry, centered away from the origin. 

Conjecture 7.4. Let f ~ 1. Given a hyperbolic equilibrium of a related averaged 
system, there is a wedge-shaped region of the form 

V = {{e,n) G : |£| < K\fi\, K > 0, i^l near } 

such that for all (e, /i) G V with s ^ 0, the center bundle equations \7.^ have a 
hyperbolic 3— frequency epicycle manifold centered near but generically not at 

the origin. 

The remark following theorem 16. 151 is likely to hold for theorems 17.31 and 17.41 

Conjecture 7.5. Let j* > 1. Generically, there is a deleted neighbourhood W-' 
of the origin such that such that for all (£,/i) G , ( |7.i?| l has a unique hyperbolic 
discrete modulated rotating wave 'Dl ^ with Ij- — spatio-temporal symmetry centered 
at the origin. 

Conjecture 7.6. Let f > 1. Given a hyperbolic equilibrium of a related aver- 
aged system, there is a deleted neighbourhood of the origin such that for all 
(e,/^) G V-' , ]7.^ has a hyperbolic 3— frequency epicycle manifold centered at 
the origin. 
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7.2. Remarks. While numerical experiments show without the shadow of a doubt 
that spirals can anchor away from a center of inhomogeneity [8, 10, 11], this writer 
would find it very satisfying to sec this result reproduced in the laboratory. From 
a resolutely profane perspective, the Belousov-Zhabotinsky reaction appears most 
likely to yield results, as it seems the easiest to control. 

It should also be noted that the study of spirals docs not start and end in the 
plane. For instance, Comanici used the dynamical system approach in her doctoral 
thesis to study spirals on spherical domains [16]. Scroll waves, the 3-dimensional 
analogues of spiral waves, have also attracted attention from the physics and car- 
diology communities in recent years [1, 15,47,48]. 

7.3. Future and Related Work. We finish this paper with a list of problems and 
open questions: their solutions would improve our knowledge and understanding of 
spiral waves and excitable media. 

7.3.1. MRW, TW and MTW. When the FESB becomes too complex, only RW 
and TW are easily amenable to characterization. Even then, TW have not yet 
been tackled. The current approach [9-11] needs to be suitably modified so as to 
accommodate the extra variable that appears in H7.1|l and (|7.2|l . 

7.3.2. An Explicit Center Manifold Reduction Theorem. While extremely powerful, 
the CMRT has the disadvantage of being a strict existence theorem: it tells us that 
the dynamics on the center manifold are given by an ODE system with certain 
symmetries, but it does not provide the explicit relation between that system and 
the original semi-flow. 

In particular, when a specific center manifold system is studied, we have no way 
of knowing if it corresponds to a 'viable' excitable system, i.e. if it is 'attainable' 
in any way from such an excitable system via the CMRT. 

Recent observations by Lajoie and LeBlanc [33] suggest that it might be pos- 
sible to efficiently relate the coefficients of a RDS to those of the center manifold 
near a traveling wave. Is there a direct and efficient way to compute the relevant 
coefficients of the center bundle equations directly from the PDE, near any type of 
relative equilibrium or relative periodic solution? 

7.3.3. Spiral Groupings. So far, the model-independent approach based on the 
CMRT has only been used to study isolated spiral waves. Even though experi- 
ments by Li, Ouyang, Petrov and Swinney [36] have shown that spiral waves can 
be isolated with the help of a laser, they are rarely found in that state in excitable 
media (see for instance the two illustrations on p. [2). 

Spiral groupings, where two or more spirals rotate around a common center or 
one another, have much different dynamics, as can be attested by the recent numer- 
ical simulations of Pertsov and Zariski [58] . Some of the showcased interactions are 
somewhat analogous to already-obtained results about epicycle drifting, which begs 
the question: how can the current approach be altered to apply to spiral groupings 
as well? 

7.3.4. Global Spiral Dynamics. Finally, the CMRT can only be applied to local 
neighbourhoods of spiral wave solutions. Yet, many spiral interactions are global in 
nature (see [58] for details). As of now, there is little machinery short of numerical 
simulations to deal with spiral dynamics on a global level. How can this situation 
be remedied? 
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Appendix A. Selected Bibliography of Spiral Pattern Formation 

BZ Reaction and the Oregonator — [25,36,38,59]. 

Cardiac Tissue and the FHN Equations - [17,29,37,41,52,55,57]. 

Global Spiral Dynamics - [26,58]. 

Surveys - [28,53]. 

Dynamical System Approach - [4,6-11,16,21,22,34,35,44,45,56]. 
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